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PREFACE 


The tr* S* Naval Engineering Experiment Station, 

Annapolis, Maryland, is at present conducting tests on 
water lubricated bearings. High pressures, high water 

temperatures, high angular speeds, and fine clearances 
are involved in these tests. 

The authors first became interested in the problem 
of fluid flow through water lubricated bearings during 
an informal discussion with Mr, Watt V. Smith, Superin¬ 
tendent of Bearing Project, U, S, Naval Engineering 
Experiment Station, who is in charge of the tests. The 
Engineering Experiment Station has thus far largely 
confined its investigation to the problems of load 
carrying capacity and erosion effects. The fluid flow 
through these bearings appeared to be an interesting 
problem and the authors offered their services to the 
EES for an analysis of the flow, Mr, Smith made the 
full test apparatus available to the authors for their 
investigation, 

V/ith the hope of assisting the Engineering Experi¬ 
ment Station in the evaluation of test data, the authors 
chose the task of finding suitable formulae and graphical 
methods for determining the rate of flow through a 
water lubricated bearing from known pressure and tempera¬ 
ture conditions in the bearing. 

The authors wish to express their deep appreciation 
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to all who have aided in this investigation. They are 
especially grateful to Associate Professor George H. Lee, 
XJ. S. Naval Postgraduate School, for his continued 
guidance and for his careful review of the completed 
manuscript. 

The authors are also grateful to personnel of the 
IJ. S, Naval Engineering Experiment Station for their 
whole-hearted cooperation in providing the test apparatus 
and for their assistance in conducting this investigation. 
Mr. Watt Y. Smith suggested the investigation and gave the 
authors many helpful hints as the investigation proceeded. 
Eurther, an acknowledgment would he incomplete without 
mention of Mr* M. Levinsohn and Mr. W. J. Yitelozzi, who 
were tireless in their efforts to keep the hearing 
apparatus in operating condition and also to assist the 
authors in conducting their test runs. 


iii 



TABLE OF(CONTENTS 

Title Page 

Certificate of Approval i 

Preface ii 

Table of contents iv 

List of Illustrations v 

Table of Symbols and Abbreviations vi 

Chapter I Introduction 1 

Chapter II Origin of Data 6 

Chapter HI Theoretical Solutions 18 

Chapter IV Solution of EES Problem . 23 

Bibliography 45 

Appendix I Derivation of Formula for Laminar 

Flow Between Concentric Cylinders 45 

Appendix II Derivation of Formula for Change of 
Clearance With Change of Temperature 
and Pressure 50 

Appendix III Bata Sheets and Calculated Data 61 


iv 



LIST OF ILLUSTRATIONS 


Figure Title Page 

I Cross-sectional sketch of test bearing. 8 

S Photograph of test apparatus. 5 

3 Cornish’s plot of dimensionless parameters 17 

4 Page’s plot of critical angular velocity vs. 

Reynolds number. 22 

5 End view of bearing. 25 

6 Plot of kinematic viscosity of water vs. 

temperature. 32 

7 Plot of dimensionless parameters through bear¬ 

ings 1 and 4. * 34 

8 Leakage factor curve. 37 

9 Sommerfeld chart 37 

10 Plot of absolute viscosity of water vs. temper¬ 

ature at various pressures. 39 

II Plot of dimensionless parameters for a water 

lubricated bearing. 41 

12 Perspective view of a cylinder. 45 

13 End view of concentric cylinders. 45 

14 End view of bearing. 52 

15 Enlargement of stressed element in figure 14. 52 


V 



TABLE OF SIMBOLS AND ABBREVIATIONS 


A - area 

a - journal radius 

b - bearing inner radius 

C - diametral clearance 2{b - a) 

0 - bearing shell outer radius 

D - journal diameter 

E - Modulus of elasticity in tension and compression 

EES - U. S, Naval Engineering Experiment Station 

e « eccentricity ratio (ratio of the distance between the 
centers of the journal and bearing to the radial 
clearance)• 

F - force 

g - acceleration due to gravity 

h - radial clearance (b - a) 

K - Taylor*s criterion of instability 
L - length 

M - mass rate of flow 

m - hydraulic radius 

N - angular journal speed (RPM) 

P - pressure difference between two points 
P” - bearing load per unit of projected area 

“• equivalent bearing load due to side leakage 
Px - pressure at any point x 
0, - volume rate of flow 

q - volume rate of flow due to pumping action 
R - Reynolds number 

r - radius 

T - temperature difference between two points 


vi 



- temperature at any point x 

IT - average fluid velocity in azial direction 
7 - displacement in radial direction 

W - displacement in axial direction 
Z - axis parallel to shaft 

c^UjZhy coefficient of linear thermal expansions of shaft, 
bearing metal and bearing backing, respectively 

P - constant in Suzuki^s equation for turbulent flow 

y - Poissons ratio 

AC ) "■ small finite difference of quantity in parantheses 
S - unit increase in volume 

- unit deformation in radial direction 

6t “ unit deformation in tangential direction 

- unit deformation in axial direction 
yW - absolute viscosity 

yy - kinematic viscosity 

- density 

dp - normal stress in the radial direction 
Ot - normal stress in the tangential direction 

- normal stress in the axial direction 
7r - shear stress 

(P - angle 

- angular Journal speed 

” critical angular velocity 


vii 



CHAPTER I 


INTRODUCTION 

Water lubricated bearings are used in a number of 
applications, the most important in Naval Engineering 
being the stern tube and strut bearings for ship pro¬ 
pulsion machinery. These bearings, as well as the 
majority of water lubricated bearings, employ large 
clearances and comparatively slow journal speeds. 

The analysis of the above mentioned bearing is quite 
different from the problem posed by a full journal bearing 
with small clearance, high shaft speed and with lubrication 
by water at high temperature and pressure. Such a bearing 
might find application in modern high-speed steam turbines 
and pumps where the use of water as the bearing lubricant 
would simplify the present elaborate gland seal system 
and make iiiq)OSsible the introduction of oil into the 
feed system from turbine or pump bearings. The analysis 
of axial flow of water through such a bearing, subject to 
light loads, is the subject of this investigation. 

Analyses of ’’end leakage” or axial flow, in oil 
lubricated bearings are given in various lubrication 
publications. These analyses, as far as the authors have 
been able to determine, confine their treatment to stable 
or laminar flow, and are therefore inadequate for our 
problem where large pressure drops through the bearing 
and high speeds of journal rotation will necessarily 
result in turbulent flow. 
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Extensive investigations have been made of water flow 
between rotating concentric cylinders* Neglecting eccen¬ 
tricity, our bearing problem reduces to the particular case 
of water flow between two concentric cylinders, the inner 
cylinder rotating and the outer cylinder fixed. While 
an analysis of this case has been made by a number of in¬ 
vestigators, the tests at the EES are with smaller clear¬ 
ances and higher speeds of journal rotation than any re¬ 
sults known to the authors. 

The problem of stable flow between rotating con¬ 
centric cylinders is quite involved and there is consid¬ 
erable disagreement among investigators as to what cri¬ 
terion to use to determine stability. The most widely 
accepted criterion for stability of flow between rotating 
concentric cylinders with no axial flow was developed 
mathematically and verified experimentally by Taylor (18)* 
Since Taylor first published his criterion in 19£3, a 
number of investigators have verified his results. Eor 
the specific case of only the inner cylinder or journal 
rotating, Taylor finds stable flow existing below a 
definite and determinable speed of rotation, while above 
this speed unstable or turbulent flow can be expected* 
Goldstein (s) verified Taylor’s results and further 
found that the same criterion was satisfactory for axial 
flow with low values of Reynolds Number, R. Eage (?) 
extended Taylor’s and Goldstein’s investigations for 

^Numbers in brackets refer to bibliography 
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values of R up to 680, Fage*s results are applicable to 
our bearing problem and are reproduced in figure 4, 

The criterion evolved by Taylor and extended by 
Goldstein and Fage appears satisfactory for determining 
the conditions which must be satisfied for stable flow 
in a water lubricated bearing; but, of course, it furn¬ 
ishes no solution for the rate of flow through the bear¬ 
ing. Empirical formulas as well as graphical solutions 
based upon dynamical similarity exist for turbulent flow 
between concentric cylinders with the inner cylinder 
rotating. By taking into account the effect of eccen¬ 
tricity caused by the load on the bearing, we may apply 
the above mentioned formulas and graphical solutions to 
our bearing problem. 

In the following chapters the authors will present 
both empirical equations and graphical solutions for 
flow through water lubricated bearings. In this in¬ 
vestigation, however, the authors find that the solu¬ 
tion is adequately represented by graphical plots such 
as shown in figures 3 and 11. Furthermore, graphical 
solutions may be used for both turbulent and laminar 
flow. 

In analyzing a problem of this type it may be 
assumed that at sufficiently slow speeds of journal ro¬ 
tation and with moderate pressure drops through the 
bearing, laminar flow will result. In the classical 
equation for laminar flpw between concentric cylinders. 
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the rate of flow varies directly as the clearance cubed 
and the pressure gradient axially along the cylinder, 
and inversely as the viscosity of the flowing fluid. 
Therefore, accurate knowledge of the running clearance 
of the bearing is essential for predicting the flow. 

The temperature conditions throughout the bearing must 
be known because of their effect upon the viscosity 
of the water and, to a lesser degree, upon the clearance 
in the bearing. Accurate knowledge of entrance and 
exit pressures in the bearing is also essential. These 
pressures affect the bearing clearance and, most im¬ 
portant of all, affect the rate of flow. 

The EES bearing poses a number of problems which 
have not ordinarily been present in previous experiments 
to determine the flow between rotating concentric cyl¬ 
inders. 

The first and most difficult problem is accurate 
determination of the running clearance of the bearing. 

As the mean diametrical clearance as measured at 
68° E, is only ,0038 inches, any small variation in 
this clearance results in relatively large change of 
flow. There is a difference in coefficients of thermal 
expansion of the journal and the materials composing 
the bearing, which results in a change of clearance 
with change in temperature conditions in the bearing. 

At high temperatures there will necessarily be a temp¬ 
erature gradient across the bearing which further 



complicates the determination of the clearance. Further, 
in order to withstand the high pressures involved and 
at the same time prevent leakage, the design and assembly 
of the EES bearing required certain end restraints upon 
the individual bearings in order to retain them in place. 
The magnitude of these end restraints are unknown and 
their effect upon the clearance cannot be determined. 

Then too, there is a light load on the EES bearing 
which makes necessary the consideration of the effects 
of eccentricity. As the eccentricity of any particular 
bearing will vary with changes of load, Journal speed, 
and pressure drop across the bearing, the exact solution 
of this problem is extremely difficult. 

Under conditions of low rates of flow and high 
Journal speeds, the temperature of the water increases 
as it flows through the bearing. The resulting change 
in viscosity may have an appreciable effect upon the 
flow through the bearing. 

The problems posed by the EES bearing will exist 
in varying degrees in any operating water lubricated 
bearing. Detailed analysis of the EES bearing is made 
in Chapter IV, 
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CHAPTER II 


ORIGIN OF DATA. 

The test hearing consists of four hearings on a 
vertical test shaft of constant diameter as shown in 
figure 1* The shaft is necked down to 3/8” diameter 
before extending through the upper and lower stuffing 
boxes in the hearing housing. The two central hearings 
are in a common hearing shell. All hearings have 1/8” 
radius fillets on each end of the hearing metal. The 
central hearings can he loaded in a radial direction 
by means of a spring. The shaft is rotated by means 
of a belt driven pulley; the pulley being mounted on 
the shaft above the hearing. Combination sleeve and 
thrust hearings are installed on either side of the 
drive pulley to carry the weight of the shaft and to 
prevent any radial thrust from being transmitted to the 
test hearing. For convenience the hearings have been 
numbered from one to four from top to bottom. The 

i 

bearing housing is supported at its base on a rigid 
stand, A photograph of the test apparatus is shown in 
figure Bo 

The finish of the journal varies from 2,0 to 3.25 
micro inches (RMS) and the finish of the bearing metal 
surfaces varies from 32 to 50 micro inches (RMS), 

Each of the four bearings may have different metal 
temperatures. The temperatures of bearings 1 and 4 are 
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recorded as and tg, the temperature of the middle two 
bearings as tig. OSiese temperatures are measured by- 
thermocouples. Temperatures tg and are taken at the 
top and bottom of -the bearing inner casing and tempera¬ 
tures til, tig, and tiQ are taken on the outside surfaces 
of the bearing housing near the top, middle, and bottom 
respectively. All thermocouple positions are indicated 
in figure 1. 

All temperature measurements are taken by means of 
iron-constantan wire thermocouples and are recorded by 
a continuous-reading recorder. Ambient temperature is 
recorded as ti 5 . The temperatures obtained by these 
thermocouples are correct to plus or minus 2 deg. i*. 
in absolute temperature and plus or minus 1 deg. F. for 
temperature differences. 

All pressures are determined by means of Bourdon 
tube pressure gauges. These gauges are connected to the 
water inlet pipe, both discharge pipes, and to the space 
behind bearings 2 and 3. The gauge locations and pressure 
ranges are also indicated in figure 1. The error of-ftie 
pressure gauge readings does not exceed 1^^ of the maxi¬ 
mum reading to which the scale is graduated; except that 
the error in reading the graduations between 33-1/35^ and 
66-2/3" of the maximum scale reading does not exceed 1/2^ 
of the maximtim scale reading to which the scale is grad¬ 
uated. 
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The mean diametral bearing clearances for each of 
the four bearings are designated as Cg, Cg, and 
respectively. All reference measurements are taken un¬ 
der standard conditions at ee^F, The accuracy of the 
measurement of the journal and bearing diameters is 
correct within plus or minus 2 ten thousandths of an 
inch. 

Water is admitted through a pipe to the circular 
groove between bearings 2 and 3 at pressure p^^ and 
temperature tg, flows axially through bearings 2 and 3 
to the spaces between bearings 1 & 2 and 3 & 4* These 
spaces are interconnected by approximately 1 / 8 ” clear¬ 
ance in back of the bearing shell which houses bearings 
2 and 3, so that the pressures will be the same, pgj 
however, the temperatures need not necessarily be the 
same and we use tg and tg to represent the temperatures 
in the upper and lower spaces respectively. The flow 
through bearings 1 and 4 is again axial and the result¬ 
ing pressures and temperatures after flow through the 
bearings are P 4 , t^ and p^, t 3 at the top and bottom 
respectively. The discharge from each end is led to 
separate collecting tanks located on platform scales. 
Here the quantities of flow are weighed and tempera¬ 
tures are taken as t^^^ and t^g for flow from top and 
bottom bearings respectively. To insure that the dis¬ 
charge pressures pg and P 4 are the same a mercury man¬ 
ometer is connected across the discharge lines. jFlow 
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from either line can be throttled by means of a globe 
Talve to equalize the merctiry col\imns of the manometer. 
The materials of construction with approximate per¬ 
centage compositions are shown below: 

Shaft TJ.S.I.S, Symbol 440c alloy steel 

(loO C, 17 Or, 0o75 Mo) 

Bearing Material- Beariiua #10 (70 Cu, 20 Pb, 10 Sn) 

Bearing Shell- Fed, Specs, 46S1-M Steel 

((6.310,0,75 Mn, 0,25 Si, 0.35 Cu, 

0,25 Ni) 

Inner Casing- COMP M Bronze (88 Cu, 6,5 Sn, 4 Zn, 

1,5 Pb) 

Outer B;ouslng- Fed,Specs, 4952-AN Forged Steel 

(0,45 C, 2.75 Ni) 
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CHAPTER HI 


THEORETICia. SOIXJTEONS 

In this chapter the authors vd.ll present general 
solutions for the rate of axial flow through a water 
lubricated bearing. Their solutions will include both 
laminar and turbulent flow and they will also discuss 
the various criteria for stability of flow. 

As previously stated in Chapter i, laminar flow 
may be expected in a water lubricated bearing under cer¬ 
tain conditions. The equation for laminar flow between 
two concentric cylinders is given by Lamb (is) as 



Where 

Q is the volume rate of flow; 

P is the axial presstire difference through the bearing; 
^ is the absolute viscosity; 

L is the length of the bearing; 
b is the inner radius of the bearing; 
a is the radius of the journal. 

Equation (l) may be put into a more convenient form 
by substituting for b its equivalent, a f h, and expanding 
to* 



* See Anpendix (I) for mathematiceLl derivation. 
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1/Vhen the radial clearance is small compared to the 
radius, equation (l) hecomes 

Q. = 

6^ L (3) 

Equation (s) above is also the fonnula for viscous 
flow between two parallel flat plates. 

The diametral clearance 0 is frequently used in 
bearing calculations. Substituting C for Sh and D for 
2r, equation (3) becomes 


96 jU z, 



Further, the eccentricity of the journal must be con¬ 
sidered, With the journal stationary and for laminar flow 
the eccentricity may increase the quantity of flow for a 
given pressure difference by as much as 150 percent over 
that with no eccentricity. The effect of eccentricity 
may be obtained by reference to Orloff (14)• Equation 
(4) then becomes 

9d ■ (1+ 1-5*2) (5) 

where ® is the eccentricity ratio and has a maximum vdlue 
of unity when the journal is hard over against the bearing. 
The mean axial velocity U may be obtained from equa¬ 
tion (5) above by dividing Q, the axial flow, by the cross- 

sectional area of the flow path, j’ormula (5) 

2 

*Approximate area, where radial clearance is small compared 
to the journal radius. 
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then hecomes 


U= ^ — (l-^l,5e2). (6) 

48yC/ L 

Equation (5) does not consider the effect of rotation 
-with an eccentric journal. Orloff gives a formula for 
additional end leakage due to the action of pressure de¬ 
veloped in the loaded zone. With approximations this 
formula reduces to 


1 = 


X 10ND® (C/D)^*^ 

D/L(l i- 0.6D/L 
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Where 

q is volume rate of flow (ft®/sec); 

N is angular journal speed (RPM); 

D is journal diameter (inches); 

is angular journal ^eed (rad/sec); 

P”is projected hearing pressure. 

An examimtion of equations (5) and (7) will show 
that for light loads and reasonably large pressure drops 
through the bearing the additional flow as computed by 
equation (7) becomes negligible small. Therefore, equa¬ 
tion ( 5 ) is sufficient for effecting a solution of laminar 
flow in the EES bearing. 

Thus far entry and exit losses have been neglected. 
Where the pressure drop used in effecting a solution of 
the flow problem is obtained from the difference of 
pressures measured in the high and low pressure chambers, 
at the entrance and exit of the bearing respectively, the 
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effect of these losses is to decrease the measxired 
pressure drop. Cornish (s), in his experiments with 
flow between concentric cylinders with the inner cylinder 
rotating, found that the speed of rotation had very little 
effect upon the entrance and exit losses and that these 
losses approximated Cornish attributed his low value 

2g 

of these losses to slight ’’bell mouthing” of the entrance. 

The bearings used by the authors were rounded at the 
entrances and exits by a 1/8 inch fillet. Their apparatus 
provided no means for determining entry and exit losses. 
The authors found, however, that the actual rates of flow 
in the laminar flow region were generally as great or 
greater than the calculated flow as determined by neglect¬ 
ing these losses. Furthermore, assuming these losses 

are approximately u£ , then a correction of mav be 

2g 2g ' 

made to the measured pressure drop across the bearing. 

This correction becomes appreciable only at very high 

velocities. 

No such simple formulas as previously given for 
laminar flow will suffice for turbulent flow in a water 
lubricated bearing. As is usually done in turbulent flow 
problems, the principle of dynamical similarity may be 
applied. If we follow the methods used by Cornish we 
obtain these dimensionless parameters and their relations: 


Pm^ 

JJ.'m 


* ^ * 

m 


'SiCO 

U > 
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m UPq 
_ 


mU 

7^ 


m 



( 9 ) 


where 

m is the hydraulic radius; 

^ is the kineiuatic viscosity; 
r is the mean radius; 

UOq is the critical angular velocity. 

Cornish was able to reduce the number of parameters 
by multiplication of individual parameters, thus 


mU . r • mCtP * 

7^ m U 


( 10 ) 


/rm — - 


( 11 ) 


where is the density of the fluid, 

Cornish then plotted m^ against R for various 

values of vmUfi to obtain his figure 7, as reproduced in 
figure 3. Figure 3 is a very useful family of curves and 
provides a means for quickly obtaining a solution for the 
flow through a water lubricated bearing. This solution 
neglects the effects of eccentricity. Cornish used values 
of r/m ranging from 341 to 624 to obtain his curves and 
states that valid results can be expected over this range. 
A large number of bearing problems might be expected to 
fall within this range. 

A brief mention should be made of Cornishes apparatus. 
He experimented with five sets of apparatus, the general 
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construction "being a smooth brass cylinder rotating in¬ 
side a hollow cylinder also of brass and finished smooth. 
Rotation was accomplished by means of a D.C, motor, giv¬ 
ing a range of speeds from 0 to 2000 RRM, The mean radius 
r varied from 6 to 10 centimeters and the length of cyl¬ 
inders varied from 15 to 28 centimeters* His values of- 
varied from 138 to 624* 

The values of r/m used by the authors were much 
higher than those used by Cornish, varying from 1319 to 
1369* However, Cornishes work appears more closely 
allied with the investigation of the authors than the 
work of other experimenters* It is of considerable in¬ 
terest to note Cornish’s curves converge rapidly in the 
highly turbulent region and are thus essentially inde¬ 
pendent of shaft rotation in this eo’ea* This is import 
ant in this investigation for the hi# journal speeds 
and large pressure drops involved necessarily result in 
highly turbulent flow* 

It should be noted here that the results obtained 
by Cornish do not entirely agree with results of later 
investigators, namely Goldstein (8) and Fage (?)* 

Cornish differs in his methods of determining critical 
angular velocity for the breakdown of stable flow* 

The criteria for stability used by the various investi¬ 
gators will be discussed shortly; but regardless of his 
determination of the beginning of instability, the major 
portion of Cornish’s curves* should be reliable over the 
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range of r/m from S41 to 624. ^ile the authors used 
values of r/m varying from 1319 to 1369, they obtained 
fair agreement with Cornish as is shown by comparison 
of figures 3 and 11, 

It must be eu^hasized that figure 3 neglects the 
effects of eccentricity and that eccentricity is neg¬ 
ligible in figure 11 except the curve for "0, 

for this curve with the rotor stationary the apparent 
eccentricity ratio as calculated by formula (5) varied 
approximately from 0.7 at R equal 100 to 0,3 at R equal 
320. With light loads on the bearing and at high journal 
speeds the eccentricity is negligibly small and figures 
3 and 11 will give adequate solutions to applicable 
bearing problems. For values of r/m near 500 Cornish’s 
plot as shown in figure 3 is recommended for bearing 
problems. With values of r/m near 1345, the plot 
shown in figure 11 may be used, ^ 

A number of formulas for turbulent flow exist which 
are applicable to the flow through a water lubricated 
bearing, Davies and White (4) give a general formula 
for turbulent flow in a pipe of euay form with smooth 
boundaries, which when applied to this problem gives 

Ssii * 0.0283 ( ^(12) 

yMin \ ^ J 

Suzuki (16) derived his formula for flow between 
concentric cylinders, the inner cylinder rotating, by 
working from Karman’s law of velocity distribution. 


19 



For our problem his formula reduces to 



Suzuki did not give a value for the constanty5 which 
appears to vary with clearance, Cornish found close 
agreement with/^ equal to 0.8 for values of R from 
1000 upward. The authors, in their investigation, 
found insufficient agreement with Suzuki*s formula to 
merit quoting any value of the constant ^ . 

The problem of stability of flow of a viscous fluid 
between two infinitely long coaxial rotating cylinders 
with no axial flow was first solved mathematically and 
verified experimentally by Taylor. Taylor*s formula for 
our case of the inner cylinder only rotating is 

_(aJl -kL - K (14) 

i? 


Where K = 0.0571 ( 1- 0.652 ^-^0.00056 1 - 0.652 • 

Equation (14) above was rewritten by Fage in the more 
convenient form 



(15) 


Therefore, for the case of no flow; i.e., R equal to zero, 
the speed of rotation where turbulence first sets in is 


20 



easily calculated by using formula (15) above* 

Goldstein investigated the stability of flow under 
the influence of a pressure gradient parallel to the axis 
and with the inner cylinder rotating at a uniform speed. 

He found Taylor*s formula gave good results at low values 
of R; i.e., R < 10. Fage experimentally extended Gold¬ 
stein’s work, obtaining a plot of 2 (flgaTn(a -tb) versus R 

m , 

for values of R up to 680. Taylors point for R equal to 
zero fits smoothly on Fage’s curve. Fage’s plot is shown 
in figure 4. 

Cornish found essentially independent of 

r/m over a range of values of r/m from 341 to 6S4. An 
examination of equations (14) and (l5) shows that accord¬ 
ing to Taylor’s theory for no axial flow r ^c varies 
as (r/m)^. Since Goldstein and Fage later found that 
Taylors criterion for stability might be modified to'^xial 
flow considerable doubt has been cast upon the method 
used by Cornish to obtain the critical rotational ^eed. 

Fage’s plot appears most useful for determining the 
breakdown of stability for flow through a water lubricated 
bearing, but it offers no means of determining the rate 
of flow. The authors feel that a graphical plot such 
as made by Cornish represents the best solution to our 
problem of flow through a water lubricated bearing* 
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CHAPTER 17 

SOECJTI ON OE OHE EES PROBLEM 

The following particulars are given for the EES 
problem. 

Measurements (68°F.) 


Main bearing, length 

4,75 

in. 

Central 

circumferentialgroove, width 

.75 

in. 

End bearing, length 

2.00 

in. 

Journal, 

diameter (mean) 

2.5003 

in. 

Bearing, 

mean inside diameter, No.l 

2,5041 

in. 

Bearing, 

mean inside diameter, No.2 

2,5039 

in. 

Bearing, 

mean inside diameter. No.3 

2,5044 

in. 

Bearing, 

mean inside diameter. No.4 

2.5041 

in. 

Bearing metal, thickness 

.250 

in. 

Bearing, 

mean outside diameter, Nos, 1 & 4 

5.2493 

in. 

Bearing, 

mean outside diameter, Nos. 2 & 3 

4.8750 

in. 


Materials of construction 

See Chapter II, Origin of Data, 
Physical nronerties of materials 

oC 

(in,/in./°I',) 

Journal 5,9 x 10"® 

Bearing metal 10.28 x 10"® 

Bearing backing 6.51 x 10"® 

Load applied to bearings 
Load applied to each bearing 
Temperature, water 
Pressure drop across end bearings 


(psi) 


Compressive 

(r) 


10 X 10® 12,000 

28.8 X 10^ - 

100 lb, 
50 lb, 
83-20EOr. 
49-1120 psi 


.33 

.26 
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Angular journal speeds 0 - 12,000 R.P.M. 

Calculations are made for "bearings 1 and 4 only. 

As shown in figure 1 there is a radial clearance of 1/8 
inch behind bearings 2 and 3 permitting exchange of flow 
between the discharges of these central bearings. Due to 
this construction, the actual measurement of flow through 
these individual bearings is impossible. This design 
feature provides equal pressures at the entrances of bear¬ 
ings 1 and 4. 

For laminar flow through the EES bearing the quanti¬ 
ty of flow may be readily obtained by application of e- 
quation(5). When the flow is turbulent we must either 
obtain a solution from a graphical plot of dimensionless 
parameters or apply empirical equations of Icnown accuracy. 
Inasmuch as a graphical plot of dimensionless parameters 
will give a solution for either laminar or turbulent flow, 
the authors believe the graphical method offers the best 
solution to the EES problem. 

Sufficient data was taken and a plot of - 

versus R for constant values of rmOi^ is shown in 
Figure 11. It should be emphasized here that Figtire 11 

is theoretically applicable only to a bearing with an ^ 

m 

approximately the same as the EES bearing, i.e., 1345. 

Similarity would further require the same eccentricity 

ratio. Ooncerning ^ , the authors agree with Cornish 

that exact equality of £. for similarity is not required 

It 

and therefore figure 11 may be used for a fairly wide 

range of values of , Eor values of £ near 500 

m m 
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Cornisli*s curves as shown in figure 3 are recommendedo 

The solution of the EES problem may be best illus¬ 
trated by calculating one of the points in figure llo A 
calculation of run 6* will illustrate the method followed 
by the authors in obtaining points for the curves of 
figure 11, The foot-pound-second system will be used 
in the following calculations unless indicated otherwise. 

The factors which must be considered in determining 
the mean clearance are: (l) the expansion of the bearing 
and shaft due to change in temperature from 68°F; (e) 

the pressure effects on the bearing surfaces; and (3) 
the increase in shaft radius due to rotation. 



Figure 5 End View of Bearing 

Considering first the expand on of the bearing due 
to temperature we wish to compute the displacement of the 
inner radius b. V7e may consider the bearing to be made 
up of two shells shown as parts 2 and 3 in figure 5 above. 
Part 2 is bonded to part 3, Upon heating the bearing 
there will be created a pressure between the bearing 

*See Data Sheet and Calculated Data 
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metal and the hacking since e^p is larger thanoC , The 
radial stresses in parts S and 3 due to the difference in 
coefficients of expansion are 

(r% a 

and 

(Tz = 

where p is the pressure at the bonded surface. 

The pressure p is determined from the condition that 
during heating both tubes have the same circumferential 
elongation, i.e. 



Substituting for p in the expression for (r~- gives 



Substituting numerical values in the above expression 
we have 

-828 psi 
compressive. 


(TZ = (10.gS - 6.51)10"° (92 - 68)10 x 10^ . 

1 I .25 / lO X 10^ \ 

• 94 V28,8 X 10®/ 
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Some stresses may exist in the parts of the bearing 
due to the bonding process* We have no means of measur¬ 
ing these stresses which are normally tensile stresses 
in the bearing metal. As the stress due to pressure as 
computed above is only 828 psi compressive and as the 

bonding stresses are normally tensile, it is to be ex- 

/than 

pected that the combined stresses will be less'the yield 
strength of the bearing metal, i.e., 12,000 psi. There¬ 
fore, no plastic flow of the bearing metal is to be ex¬ 
pected. 

The unit deformation in the tangential direction is 


^2t r - 68) 

^2 


Therefore 


Ad 


Acircumference 
2 W 


2 CTd a 


21 


z7r 




Substituting numerical values gives 


Ad 


10.28 X 10”®(92-68) - I 1*5 = 244 x lO^^in. 

I. 


The unit deformation in the radial direction is 

£gr : oCg(t - 68)- )rpt 

% 

= 10.28 X 10"®(92-68) - = 273 x 10"®in./in. 


The change in radial thickness of the bearing metal is 



Hence 


. Ad - A(d - b) 

= j£44 - 683 10"® z 176 x 10"® Inches< 
which is the change in b due to change in temperature» 

The effect of the pressure on the bearing shell will 
be considered next. Due to the construction of bearings 1 
and 4, the pressure pg will be present on the outer surface 
of the bearing shell since the gaskets to prevent water 
leakage in back of the bearings are on the down stream 
ends of the bearingso 

Formulas (28) and (52) of Appendix II will now be 
employed to determine the radial displacement, Ab, of 
the bearing inner radius. 

The actual bearing is composed of bearing metal ,25 
in, thick with an S of 10 x 10® psi and the bearing back¬ 
ing is ,94 in, thick with an 1 of 28,8 x 10® psi so that 
negligible error will be introduced by assuming a homo¬ 
geneous steel bearing for calculating Ab due to pressure 
effects only. 

Since the bearing temperatures t^ and t 5 and the 
outer housing temperatures tj^Q, and t^^g correspond, 

the temperature gradient aspects of formulas (28) and (32) 
will be omitted. Equation (52) then reduces to 
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where ^end pressure, pg, acting on each end of the 

bearing and p^ may be taken as the arithmetic mean of pg 

and P 4 * Then upon substitution we have 

255 _ 2 X .26 ({1.25)^ 153 -(2.44)2 eSs ) 

" 28.8 I 10® 28.8 X 10® ((2.44)2 - (l,g5)2j 

= -3o6 X 10”® in./in. 


Substituting the above value of 6^ in equation (28) of 
Appendix II g^ves 

Ah = (1+r) 


r(l -2jr) - c^PoJ 


E 


>2 M hS 


r(o2 - b2)l E 


- jy- (-3.6 X 10“®) 
i--r f 


Whicl^ upon further substitution gives 


Ab = (It.26) 


1.25(1 - 2X.26 ) /( I.25)^ x 153-(2.44)2 x 255 . 


28.8 X 10® 


(2.44)2 « (1,25)2 


(1.25)2 (2.44)2 

f 

'153..255 

1 1.25 X .26 (-rr A y 4 0-6) 

1.25 ((2.44)2-.(l.25)2j 

1 

^^28 .8 X 10^ j 

1 + .26 


Ab = -14 X 10“® in. 

which is the change in b due to pressure. 

To calculate A. a for the shaft it is assumed that the 
shaft is at a uniform temperature throughout. Therefore 


Aa = (t - 68) 

- 1.25 X 5.9 I 10”® (92-68) r 177 x 10"® in. 

The effect of centrifugal force on the shaft radius is 
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negligible for all runs of the experimental data. Assuming 
the worst possible case, i,e., the shaft to be a thin disc 
rotated at 18,000 rpm, then* 

Aa a (1 - ^ 2) /f . ^ g a? . , .. 

gH, 

® ( 1 . 85)3 

*fl “(.26)g) (*89) - 

^ 32.8 z 12 I 88.8 z 10® 

z 14.9 z 10“®ln. 

The^a calculated above due to centrifugal force assumes 
the worst possible condition used in this experiment and is 
relatively small in comparison to the tolerance for reference 
measurements; therefore, the effect of centrifugal force on 
the change of clearance will be neglected. 

The total change in radial clearance is the sum of the 
changes of /ib due to ten^jerature and pressure, and Z^a,i,e. 

Ah s [^-14 +176 -177^10“® in-, 
a -13 X 10“® in. 

The diametrhl L clearance at 68*^?, is 

C reference s 8.5041 - 2.5003 s 3800 z 10”® in. 

The actual clearance is therefore 

^ actual = [ 38 OO - 2 X I 3 J 10”® s 3774 z 10”® in. 
s 5.37 z 10”9 in.9 

See reference ( 20 ). 
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The measiared mass rate of flow through bearing 1 is 
- 25 Ibs./S min* ** -15«8 sec* * 11*01 lb •/min. 
Correspondingly, the mass rate through bearing 4 is 

“4= 1 mlnMue sec. ' 13.48 lb ./min. 

The velocities through bearin^^ 1 and 4 are* 


Un - . 144 :,, 3 C.., 2 _ . 28^7 ft. /sec. 

60 X 2*5 ^ X 3772 x 10“6 


and 

U = 13>48 X .01611 X 144 x 2 
^ 60 X 2.5 X 3772 X 10“® 


35.2 ft*/sec« 


The Reynolds Nmnbers through bearings 1 and 4 are 


R = l£-. • CTT 

R = 5772 X 10*^ X ,28.7 - 278 

1 ” 4 X 12 X 8.1 X 10”^ 


and 

R - 3772 X 10~^ X 35,2 - 

2 " 4 X 12 X 8.1 X 10~® 


The rotational flow dimensionless parameter is 


rm (jJ) 5 0, since Ct^ r 0. 

The third dimensionless parameter is 

m^ = 5.57 X 10"0- X .01611 x 204 x 52.8 

^2^L 64(8.1 X 10“S}2 X 2 

* See reference (ll) for specific volume. 

** See figure (6) for kinematic viscosity. 
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The other points for plotting the curves shown in 
figure 11 were obtained in a similar manner. It may be 
noted that no consideration of temperature gradient across 
the bearing was made in the above computations. Temperature 
measurements, shown in Data Sheet and Calculated Data, 
indicated this gradient was negligably small for all the 
test runs made in this investigation. Further, it was foxind 
that the change of temperature of the water flowing through 
the bearing was small, being a maximum of 4‘^F. for all runs 
except no. 50. The indicated rise of S^F, for this run 
appeared to be erroneous. IJ'iTith these small temperature 
rises, the authors felt justified in using the arithmetic 
mean t'lmperature of the water flowing through the bearing 
for determining the viscosities used in their computations. 

The authors were unable to obtain equal flow from 
bearings 1 and 4. The apparent trouble was that bearing 
4 was slightly skewed or restrained in the inner housing. 
This may best be shown by figure 7 and the forthcoming 
explanation. 

Figure 7 is a plot of dimensionless parameters for 

rm^ equal to 0 and 440 for bearings 1 and 4, For compari- 

son, a plot of Cornish's curve for rm^e/^ equal 0 is shown 

P' 

by a dashed line. The lower or straight portion of this 

line is the theoretical plot for laminar flow as given by 

equation 4, As the majority of the runs for rmtf'r 440 

P' 

were made at approximately 4000 RPM, we find 
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lOOjOOO 


m 

PLOT OF DiriEN5iONL£55 PARAMETERS for FLOW 
THROUGH BEARINGS 1 and 1 for And^^O 


#1 


10,000 


/ 

/ 
















2(t/am ( a+b l 17.2 

^ \yn I 

for this speed. Referring to figure 4, we find this point 
well below the c\irve and thus the effect of this speed of 
rotation on turbulence may be neglected. Therefore, if 
the bearings and .journal are concentric , we can expect the 
curves for rmu P:; 0 and rmo /a 440 to coincide. 

The curves for bearing 1 do not coincide at the 
lower pressures due to the 100 lb. load on the middle 
bearings tending to cause eccentricity and resulting in an 
increase of flow when the rotor is stactionary; while at 
4000 RPM the effect of this load on eccentricity is neg¬ 
ligibly small.* As the pressure drop and resulting flow 
is increased the eccentricity decreases as is to be ex¬ 
pected, and the curves for bearing 1 coincide. 

An inspection of figure 7 shows that the curves for 
bearing 4 are very close at the lower Reynolds Numbers but 
do not coincide at any point. Furthermore, the curve for 
vmU ^s 440 is slightly lower than the ciarve for rmU /’s 0, 
indicating an increase of flow with rotation* Such re¬ 
sults are not normal for an unrestrained bearing as the 
eccentricity should decrease with rotation resulting in 
a decrease of flow as- shown by equation (5) in Chapter III. 
The increase of flow with rotation can then be explained 
by the pumping action of an eccentric Journal as shown 

*A full discussion of eccentricity follows shortly. 
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by equation (7)* Therefore, bearing 4 must be restrained 
or skewed in its housing. 

Although extreme care was exercised in the manufacture 
and assembly of these bearings it is virtually impossible 
to insure concentricity of upper and lower bearings in 
this type of apparatus unless the bearings are concentri¬ 
cally bored in place. Due to the time required, the faults 
in bearing 4 were not corrected. It was decided to in¬ 
clude only the results obtained from bearing 1 in the 
final results of this investigation. Figure 11 is a 
dimensionless plot of the results obtained from bearing 1. 

Let us now consider the effects of eccentricity on 
our problem. By equation (5), we can expect up to 150^ 
increase in rate of flow for laminar flow with an eccen¬ 
tric rotor. In the turbulent flow region the effect of 
eccentricity is much less, Cornish finding the increase 
not more than 30^ for maximum eccentricity. 

To estimate the eccentricity of the journal we may 
refer to Kingsbury’s curve of leakage factors and to 
Sommerfeld’s chart for minimum film thickness. Sommer- 
feld’s chart gives the minimum film thickness for a 
bearing without end leakage. Our problem is to obtain 
the minimum film thickness or clearance with axial flow 
or end leakage. We may do this by applying a correction 
in the form of a leakage factor obtained from Kingsbury’s curve. 
The above curves are reproduced in figures 8 & 9 respective¬ 
ly.* 

* Figures 8 and 9 are from reference (13) 
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LEAKAGE pACTOR FOR LCWP FOR BEARING5 
CF MiNIMUK: COEFFICIENT OF FRlCTiON 



^'c }jD Z.D 3.0 ■ 

LENGTH IN DlFECTfCN T nCTTIW 

AMAL i£W^7« : 


60MMERFEID CHART FOR FULL JOURNAL BCARm 
WITH UNlFORtl VldCOSirr WITHOUT LEAM&E. 























We must first obtain tbs leakage factor from figure 

80 I'ollowing Kingsbury’s practice, our bearing may be 

assumed effectively that of a partial bearing of 120® 

arc. The length in the direction of motion then equals 

TTD/S or 2,62i^# The abscissa for entry in figure 8 then 

becomes 2.62/2, or 1.31. From figure 8 we find P_is 

Poo 

about 0.33* This means that with axial flow the bearing 

will carry only about 33^ of the load that would give 

the same minimum,clearance in a bearing without axial 

flow, i.e.,end leakage causes the journal to ride lower 

in the bearing. As the load on each individual bearing 

is 50 lb, P”is 50 , or 10 psi. Poo is then equal to 

2.5X2 

10 /.33, or 30.3 psi. The apparent bearing load, Poo , 
is used instead of P” for computing the abscissa for 
entry in figure 9. 

We can now estimate the eccentricity by use of 
figure 9. For run no. 29 we have given: 
r = 1,25 in. 

^ = 1.884 X 10*3 in, 

t = 95°F. and thus yU s 10.5®X 10*® Ib-sec/in^ * 

P”s 30.3 psi 
N = 4000 RPM 

then 

(r/h)^ H _ / 1.25 ^ ^ Y 10.5fiX 10"®X 4000 

“F’’ 1^ 0.0,01.884/ 30.3 

= 6.13 

*See Fig.10 for absolute viscosity. 
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From figure 9, we find 
1 - e r 0.84 

or / 

e = 0.16 

Then, the increase in flow is given by 1.5 e2 or about 
3.84 percent. Considering the accuracies possible in the 
determination of some of the other factors in this experi¬ 
ment this small percentage is negligible. The eccentricity 
decreases fairly rapidly with increase of speed. 

Brief mention should be made of the effect of gravity. 
The static head for bearing 1 was approximately one foot, 
or .433 psi. This figure is negligibly small in comparison 
with the pressure drops encountered in this investigation# 
It may be observed that the laminar flow line obtained 
by the authors as shown in figure 11 lies somewhat below 
the theoretical laminar flow line shown in dashed lines in 
the same figure. ViThlle, at the first glance, the dis¬ 
crepancy appears rather large, actually the error is well 
within the accuracy of meastjrements. To illustrate, 
suppose we apply our limits of accuracy of measurements, 
to a specified point on our laminar flow line and see if 
it will raise the point up to the theoretical laminar 
flow line. For run 29 we have given: 

Calculated clearance, C 3768 X 10“® in. 

Limits of accuracy of meas. + 200 X 10“® in. 

Maximum possible clearance, C^^ax. 4168 X 10“® in. 
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FIGURE 11 
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C® max.-7.S4 X 10*9 in.® 


The mean velocity through the bearing then becomes 


U = 12»8 X 0.01611 X 144 X 2 

60 X 2,5 X 77- X 4168 X 10*6 


r 28,8 ft,/seCo 


The Reynolds Number, R, is then 


m N 
R = ^ 


4168 X 10“6 X 28.8 r 319 
4 X 12 X 7,85 X 10-6 


as before. 

The pressure parameter is then 


= 7,24 X 10**9 X 0,0161LX 230 X 32.2 . 1096 

64 X 61.6 X l0“O i 8 


versus the previous value of 810, Applying this value 
^ brings this point on our curve for rmca/ r 440 

well above the theoretical laminar flow line. 

Recognizing the fact that our curve as shown in 
figure 11 is slightly below the theoretical plot, the 
authors feel that the error is sufficiently accounted for 
by the limits of accuracy of the reference measurements. 
Due consideration of this deviation from the theoretical 
laminar flow line is necessary when figure 11 is used 
for obtaining the flow in similar bearings. 
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APPENDIX I 

I, - 

DERrSTATION OF FOMLA FOR LAMINAR FLOW BETWEEN 
CONCENTRIC CYLINDERS 




Referring to figures 12 and 13 above, assume steady 
flov? of the liquid straight through the annular uniform 
cross-sectional area. Taking the Z axis coincident with 
the axis of the tube, and assuming the velocity is parallel 
to the Z axis and a function of the radial distance, r, 
from this axis, the shear stress across a plane perpen¬ 
dicular to r will be: 

shear stress = abs. viscosity i velocity gradient, or 
in differential notation 




liL 

dr 


Hence, considering a cylindrical shell of fluid whose 
bounding radii are r and (r+dj.) the differential re¬ 
tarding tangential force on the two curved surfaces is 


Since steady flow was assumed, this force must be bal¬ 
anced by normal pressures on the plane ends of the shell. 
The velocity TJ is independent of the distance along 
the flow path, i.e,. 



The differential force on the cylindrical shell due to 
the normal pressures is 

“(pi • P2) 2 7rr dr . (17) 


Equating the two differential forces, we have 




3 n 2 7r rt jcir = (pi “ P2^2 TT r dr 


hence 


^ ftiZ) = (Pl - Pg) r 
V Jr/ --JZ - ■ 

The integral of equation (18) is 
<Pl -P2> ,2 

U = - --- r ^ A log r + B. 

4^L ^ 


(18) 


(19) 


Assuming there is no slippage of the fluid at the outer 
surface of the shaft or inner surface of the outer shell, 
the boundary conditions are: 

U = 0 at r ■ a 

and U - 0 at r s b . 

Substituting boundary conditions in (19), we have 
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0 =- 


^2 A In a +B 




(a) 


and 






-t- A In b + B. 


(b) 


Then, by subtracting (b) from (a), we have 

Pi -Pa 


A (in a- In b) r 




(a2-b2) 


or 


A 


- 51 -52 


in ^1) . 

Substituting the above value of A in equation (a), we obtain 
R - 5i - 5a ^3 _ pt - Ss („2 .hai _ 

Substituting values of the constants A and B in equation 
{19), we have 

P-. -P. 


U = -. r2 Pi "P2 - “ In r+-- 

4^ L ^ 4 ^ L [in ja^ 4^L 


a2-b2l, _ -^1 '■•^a ^2^ 


Pj^-Pg (a^-b2) 


In a 


which upon simplification yields 
U 


= ^1“^2 1 

(a2-r2) . ai^-b^ 

ML 1 

1 ln/a\ 

VbJ 
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In 


(f) 


( 20 ) 



Equation (20) is an expression for the axial velocity 
at any radius r. 

The volume rate of flow, Q,, is therefore 


rb 


Q = 


r’h 


TJ dA = 
a a 


TJ 2 7r T 6ir 


which upon substitution for U becomes 

Q = 


fb 

— ^ 

Pl’P2 

(a2-r2) ^ g2.b2 /lb £^1 

ln/a'\ ' &y 



ZTTv dr 


and upon integration yields 

TrtPi-Pa ) 

8 




(b4.a4) . (bg-aS) 2 

in/b) 

(.ay 


( 21 ) 


This formula is difficult to use as the numerical values 
must be carried out 12 to 15 significant figures for 
accurate results. We may simplify equation (21) by 
substituting for b its equivalent, a-f-h. Then 

^ L ' in 


7r(pi- P2) 

‘ 8 y6{L 


(a+h)2 -a^ 


(a + h)2+a2 - ia,.t..bji...«a2 

et|). 


Since 


In (1+1) = ^ - i (a]^t I (a)^ ” ^ (a]^“' 
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Then 


ft = 


h^-a^ 4 Sah-fh^ 


, (2ah-fh2) 




t - mA Id)" - i( 


- 71^(P1-P2) ah^ fe hi sl r^i 1_ /i?_L 7 

6/4.L 1^8 a'^aj 1^3^90 \ N'^''J 


and finally 


Where * is very small number the above equation reduces 


to: 


n: ^(Pi-P,)ahg _ 

6yUL 


(B2) 
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APPENDIX II 


DERIVATION OF FORliiNLA FOR CHANGE OF CLEARANCE WITH 
CHANGE OF TEMPERATURE AND PRESSURE 

* 

In the following derivation the wall of thfe bear¬ 
ing will be considered that of a cylinder which is non- 
uniformly heated thus setting up thermal stresses* It 
is assumed that the bearing is a homogeneous material, 
or if a lining is used it is so thin that its effect 
can be neglected. The effect of pressure on the in¬ 
ner and outer surfaces, and on the ends of the bearing 
are also taken into account. The temperature distri¬ 
bution is assumed symmetrical with respect to the axis 
and constant along the axis. Therefore the deformation 
is S3nnmetrical about the axis* 

A ring is out from the cylinder by two cross-sec¬ 
tions perpendicular to the axis and unit distance apart. 
It is assumed that the cylinder is of sufficient length 
so that the ends of the element remain plane* 

The Z axis is taken as the axis of the cylinder and 
W is the displacement in the Z direction. The sub¬ 
scripts r and t denote the radial and circvunferential 
directions respectively* The symbol V denotes dis¬ 
placement in the radial direction. 

The unit deformations in the three perpendicular 
directions of an element in the above ring are: 
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s constant 





The strain can be represented as functions of the 
stresses 0^ ^ (Jf ^ the temperature • Then 

^r 5 gr - £( (ri+ (/t )-^ cCT (h) 

E E \ ' 

^t = ^ _ If (Tz + (rr) -f OCT 
E E V ^ 


Where c?c is the linear coefficient of thermal expansion and 
T the temperature difference through the bearing which 
varies with the radial distance, r. 

The unit increase in volume is 

Cz ^r ^t s B (^■i- (Tt + o<^T* 

Erom the foregoing equations, the stresses are 


0~z 

Or 

(ft 


s _E_ 

1+r 

s E 

l+iT 


E 

1+r 


(6z + 

( t 
(^t-^ 


-Jl. 
1 -Er 



oCTE 

1 -sr 


y C\ oCte 

1 -Er^^ iriST 


1 


I- 

-sr 


S)_ qcte 
1 -E^ 


(c) 
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Figure 15 


Referring to figures 14 and 15 above: 

(Tt “ normal hoop stress acting on the sides and 

nni of the elemento 

(Tr = normal radial stress on the side mn. This 

stress varies with the radius r and changes by 

an amount ( )dr in the distance dr* 

dr 

The normal radial stress on the side is 

(T ^ — dr, 

''r ^ dr 

Summing up the forces on the element in the direction 
of the bisector of the angle 6.0 yields the following 
equation of equilibrium, neglecting the mass of the 
element* 

orp rdj^ -t dpdr d0 - (r + dr) d0 * 0. 
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Neglecting small quantities 

<57- r d dr . 
dr 

Substituting the stress and unit deformation equa¬ 
tions in the abovre equation of equilibrium yields 

d% ^ 1 = 1 -LI.oo ^ 

43*2 r dr pE 1 - 'jf <^3: • 

This differential equation determines the radial dis¬ 
placement 7 for any particular distribution of temper¬ 
ature. This expression may be written in the form 

r — ^ (r7)l s oO 4-^ 

dr j_ r dr j I-’IT dr 

Integrating with respect to r, we have 

4_ (r7) - l - t -E . oO Tr-fC; r . 

dr 1 - ^ 

Integrating again, yields 

T: I t Irdr+O^r+Cg J _ (23) 

J b 

The above equation contains two constants and Cg which 
must be evaluated from the botindary conditions. Calling 
the outside pressure on the bearing Pq ^i^d the inside 
pressure P^^then 

f (Jr )r - a ’ ^ 

\ r ab » 
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The negative sign is used because normal stress is 
taken as negative for compression, A general expression 
for is obtained by substituting fr, ~ r 

and the foregoing expression for V into the stress-strain 
equation for (Tygiven in equation (c). These substi¬ 
tutions give 



1 + r 




1+ y 

irrr 



f'oC 


0 b 


Tr dr X 


1-sr i-sy 


;24) 


To evaluate constants and Og the boundary conditions 
are substituted in the above general stress equation. 
The boundary conditions are: 



Substituting, we get 


E 


n-V' 1 


Pb 


“ l+n” l-y b2 

and 
"Po * 


oG Tr dr+. 


-h.-2&+ y f. 


(25) 


1-T. 


i+rj. 


If y 1 ^ 
" 1-r c2 
0 


oC Tr dr-/Xl - 
b 


1-2 r o2 




(26) 


Subtracting, we have 


(Pi-Po) (ifT) - 

■ E- - " 



Tr dr 4- ^ 

b2 


n2 
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or 


Cg s 0 ^- 
c2-b2 


- 

ic 

■(pi.-Ppl(l-t.I.). g. 1+^ 1 

Tr dr 

E ‘ l-r c2 

b 


Az7) 


Substituting the above expression for Cg in formula (25), 
we have 

n c 


-pi = _1 

1+r 


1 -_ c2[ (pi-po) (itr) ■ 1+r i_ 

v*#" o r> I T? I ^ “v- 


l- 2 r C§b 2 


E 


1-r c2 


°^Tr dr 


v/b 


1-2^ 2 


or. 


C, -ll±mi=20. 

^ E 


b2pi-o2 po) 

c2-b2 p'li-rj 


KJ 


oC Trdr-y^^. 


Substituting the values of and into the displace¬ 
ment formula (23), gives 


T = JdJT 
r i-r 


Tr dr + B(H-n(l-2r ) /b5|ip!E^ ^ 

) ^ 

fi?f) /W4U±ii)+ 




(y ' 

n 

1 i+r 

\r(c2-^2;^ 

^ 1-r 

d 


OC Tr dr 


or 
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Equation (S8) is the general expression for the radial 
displacement at any distance r and becomes A b when evalu¬ 
ated at r = b. In order to determine the radial displace¬ 
ment by use of equation (£8) it is necessary to evaluate 
the constant strain, . This will be done from the ex¬ 
pression for (Tz appearing in equation (c) by substitut¬ 
ing m s , and using equation (£8)« V/e 

may calculate the magnitude of € 2 from the condition 
that the sum of the normal forces over the cross-section 
of the cylinder must equal the external force on the end 
of the bearing due to end pressure. 

The normal force is then 




(c2-b2) 

2 . 


(29) 
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Substituting for equation (c) gives 


or = 


E 


1+r 




-X-f 

1-2^1 


f , d7 - , 
+ dr ^ 


- ocTE 

^7 1-8 y" 

J 


Further substitution for "V and --“ from equation (23) 
gives 


(^ = E 



-f. ^ 



- ^ 

1+y 

^ i-2r 

1+)" 


3,2(i-y) 





J 


oCTr dr + 


cCTrdr- ^^c^Trdry + 


( 1-8 y ) b^i-c^o^ b^o2 ^pi-po^ _ T 6z. 
'“b^/ 


E I o2.v8 y 2‘tc2-b2} I ^ y 1 tT 


+ 


((^-^i^>t^;)\oCTrdr. 


i.l-gJL) b^l-c2po\ ^ b2c2^ ^ 

e8-b2y 


E 


V 


^2 ({ 52 ^ 138 ) I E / /+T 


oCTE 

l-2r 


Upon simplification, we have 

'c 

-SX-\ oC Tr dr + ^ ^^1-0^0 j ^ 

^ c2-b2 J 

b 

oCT 


(?^=E 


(i-r)(c2-b2) 


^ToC T . oCT .1 A . 

rr-TTTTT^TT TT=2?)''"nT?) 


Zf 


1-2 y/ 


(30) 


Substituting equation (30) in (29), we obtain 
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ocjr rT 


- T^) (c2.,,e). 


(31) 


Assuming the temperature gradient is linear, then 


T a ti (1- 

o-b 


where t^ is the inside temperature and the outer surface 
temperature is zeroo 

Evaluating the integral within the integral, we have 


c pc 

r dr - ^^i 1 

c -b I 

b J h 


(r-b)r dr a ^"^i (c^.^S) 

Z 


^i (c^«-b^) ^^i b{c^»b2) 


2(c-b 
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Integrating equation (31) gives 


E 


rcCtKca.bS) _ V^H(c3^b5) ^ VcCtj ^(c2,ia8) ^ 
2 ( 1 -y) 3(1- r) (c-i) ^ 2 ( 1 - V)(c-id ^ 


y' /b^ni-o^o) -I- t 


E 


ri-2 r 


rferT) 


s£=bi_ 

2 


1 

c-b 


! ^ b(c2-b^ ))/ - ^ ^i / c^-b^ — 

\ 3 Z /j l-2y K 2 

?s^) + 


or 


^2= - ^end X F -1+ 2 (c2 f cb t-b^ ^ b ~ ] 

E 1- y |_ 3 (pipi c-b ~ 


2 Y Pi-c2 Pq) f Y'qO ti 


Elc2-bS) 




2r )(i-ar) 


} 


(32) 


fi - Ji- / 2(q?-1^) - b^l 4- Fl - 

L c-b I 3(c2-b2) /J 1-2 vL 3(c2-b2) c 


Equation (32) above is a general expression for 6^. 
Nximerical values may be substituted in this equation and 
the value of ^ 3 j, so obtained, may then be substituted in 
equation (28) for the radial displacemento The radial 
displacement includes the effects of a temperature 
gradient across the bearing and the effects of pressure 
on the bearing surfaces© 

The total radial displacement is obtained by super¬ 
imposing the values of T-obtained by equation (28) and 
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that obtained by a uniform temperature change throughout 
the bearing. This temperature change is equal to the 
change in the temperature of the outside surface of the 
bearing. This uniform heating throughout the bearing does 
not insert any stress so that the method of superposition 
is valid. 
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AFPENDIU III 


DATA SHEET CALCULATED ' DATA 


RUN NO. 

\ 
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92 
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13 
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o 
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II 

> 
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^7 
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MTA SHEET calculated DATA 
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50 
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t, 

lo? 

98 

107 

no 
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no 
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97 
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98 
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ne 

/ 09 

96 
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no 

ne 
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96 
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107 
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ta 
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93 
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95 

106 

92 

107 


ts 
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(99 
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95 
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92 
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tt 
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Z09 
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93 
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tr 
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93 

/07 


te 

108 

181 

169 

95 
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U 
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ts 

95 

98 

f6 

19/ 

138 

98 

10/ 
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35Z. 

905 

77 

/95 

J99 \ 


i 

1 

R BOT flftC. 

— 

930 

— 

/// 

1 

ZOf 

\ /ez 



5760 

1390 

3000 

398 

/6/e 

S?2 . 

6ZO 

1 


4^3 

9ZZ 

917 

995 

1065 

/ooo ! 

/ore 








-1 

1 i 

_i 

t 








1 

[ _ 


















i 



• 




—:-1 


1 

1 










































































































/ 


MTA SHEET CALCULATED DATA 


'■'O 


RUN NO. 




4/^ 

— ■ —. ■ ■ ■ ^ 

47 

4B 

47 


Rpn 

Sooo 

8 ZOO 

8000 

9b 00 

8000 

8000 

9SOO 


A 

LOO 

szs 

(>oo 

7ZS 

IZOO 

/seo 

Z05 



1^7 

ZbO 

355 

4os 

(*70 

ezo 

/oz 


/4 = A 

3Z 

so 

SO 

SO 

SO 

50 

SO 


FLOW (Lb.) 

IS 

10 

ZS 

ZS 

ZS 

Z5 

/O 


ritiEX^ 

3-15.0 

n/o.o 

Z-00.8 

/-4€./ 

/-03.S 

0-S4.0 

b 03.C, 


TiriEffl-tw 

Z-3(,.7 

O.SJ.b 

/-33.S 

/-2i.3 

0-S4.0 

0- 450 

3-S?.S 


t, 

nc 

106 

/06 

/04 

/OB 

//O 

/3<i> 


'tz 

lib 

loe 

108 

/03 

/OB 

I/O 

/34 



ns 

/OB 

J0b> 

/03 

/or 

10$ 

130 



i/b 

lOB 

707 

/Q3 j /OB 

no 

/3Z 


ts 

Hb 

106 

107 

/04 ! /08 

//o 

/33 


tfe 

IM 

108 

/Dt 

/or 

/07 

/OB 

/30 


tr 

Hb 

10? 

107 

/04 

/OB 

//o \ 

/3S 


ts 

Hb 

to? 

/or 

/04 

/OS 

no 

133 


t? 

m 

— 

104 

/o/ 

/OS 

/Ob 

/24 


tio 

/// 

JO^ 

/04 

97 

/OS 

/07 

/ZS 


t,i 

m 

/03 

/or 

78 

105 

/Ob 

/Z4 


tfz 

/// 

103 

/or 

98 

/OS 

/Ob 

/Zb 

i 

tl3 

ns 

/OS 

/OB 

/04 

/09 

I/O 

/3/ 


t)4 

H5 

/08 

/08 

104 

/Of 

no 

/3Z 


tis 

eo 

6Z 

8z 

83 

04 

05 

8Z 


tlt> 

113 

107 

/OS 

/O/ 

/OS 

/07 

/Z8 ' 











fi 

ns 

z/o 

305 

3SS 

bZO 

770 

SZ 



l/S.S 

108.5 

/Oi.S 

J03.S 

1075 

109 

^3Z.5 


R TOP 8R6 

IS/ 

Z^7 

37Z 

374 

708 

047 

6>Z 


R BOT BPt 

188 

— 

^77 

— 

835 

lO/Z 

— 



bfZ 

900 

/J80 

/4b0 

Z770 

3470 

336 



/o?o 

977 

788 

959 

967 

995 

/4?0 


















\ 











L 







1 


1 

f-- 






1 

h - 


_S_ 

^7 ^■11 























































































MTA SHEET CALCULATED DATA 


RUN NO. 

SO 

SI 

sz 

S3 

S9 

ss 

S6> 


RPn 

‘fdoo 

IOj(i,DO 

IZ, 000 

/OjBOO 

IZ, 000 

/z^ooo 

/Z.ooo 


A 

Z05 

700 

7SS 

nys 

IMS 

I90C 

/6Z0 



IS/ 

^00 

^so 

(oOO 

605 

730 

850 


^ = A 

SO 

S! 

so 

'so 

SO 

SO 

SO 


FLOVJ (Lb.) 

10 

ZS 

ZS 

Z5 

ZS 

ZS 

ZS 


TiriES'^^ 


Z~OSsZ 

z-o<,.z 

hZO.O 

hZZ.O 

/-OS.S 

1-00.9 


nMEW-^W 

Z'O^.O 

h3S.Z 

I-3S.0 

I-Oi.o 

ho^.s 

0-5Z.O 

0-99.6 


t, 

131 

IZI 

108 

! l(o 

106 

/OS 

/09 


i'Z 

130 

IZZ 

108 

1/6 

/OS 

/OS 

/O^ 


ta 

IZJ 

1/8 

/Oi 

HZ 

. HZ 

/DZ 

/OO 



/Z9 

IZI 

lor 

1/6 

/OS 

/OS 

/09 


ts 

130 

izi 

. /08 

116 

106 

JOS 

/09 


tt 

IZ9 

US 

/OS 

1/3 

/03 

/OZ 

/O! 


tr 

131 

IZI 

/08 

1/6 

/06 

/OS 

/09 


is 

130 

IZI 

108 

116 

/06 

/OS 

/09 


t? 

IZ3 

117 

/Oo 

119 

99 

98 

— j 


tio 

IZI 

H9 

lOZ 

ns 

/o/ 

10/ 

/OS 


tn 

1/8 

izo 

loz 

IM 

/o/ 

/o/ 

/Of 


tiz 

H? 

izo 

loz 

119 

/o/ 

JO/ 

108 


tl3 

IZ9 

IZI 

108 

116 

/06 

/OS 

/OS 


tif 

130 

IZ! 

108 

116 

/06 

/OS 

/OS 


tl5 

8Z 

BZ 

78 

83 

76 

78 

83 \ 


tit, 

1Z(> 

1/6 

/OZ 

Us 

/oo 

/OO 

?? 











A-A, 

101 . 

3^9 

^00 

sso 

sss 

730 

800 



130 

//9.S 

106 

1/9 

/09 

/03.S 

/OZ 


R TOP BR& 

135 

^IZ. 

3S3 

6/Z 

S33 

669 

7Z/ 


R SOT flflfr 

— 

— 

^6? 

— 

G76 

$33 

— 



733 

2090 

17 BO 

Z9Z0 

Z3S0 

Z990 

32 60 



I3ZS 

IS 03 

I97S 

/9SZ 

1991 

/9Z/ 

/9ZO 





\ 

































[ 


1 








68 











































































mTA SHEET CALCULATED DATA 


RUN NO. 

, , -—1 

51 

S8 

s? 

60 


(pZ 

63 


RPM 

12^000 

Hjoo 

6200 

11,100 

8600 

8zoo 

SZOO 


A 

/6Z0 

Z05 

Zoo 

ZD5 

ZOO 

9Z5 

665 


Pz 

955 

/OO 

105 

199 

Ml 

ZSB 

380 


^ = A 

SO 

5/ 

S/ 

5/ 

SO 

5/ 

SO 


FLOW (Lb.) 

zs 

/o 

s 

10 

IS 

ZO 

Z5 



O-SO.D 

s-^z.c 

1-97.5 

Z~96.Z 

Z-ss.o 

I-90.S 

I-Z9.Z 



O~^3.0 

Z-5<f.k 

I-IS.D 

Z~OS.5 

Z-/Z.O 

1-19.5 

l-IO.l 


t, 

105 

I5Z 

I96> 

199 

189 

996 

MB 


tz 


ISO 

195 

197 

IB7 

Ms 

M7 


ta 

101 

IH 

195 

199 

IBf 

195 

Me 



105 

196 

199 

MS 

/as 

MS 

197 


is 

}0^ 

ISO 

199 

196 

IB7 

196 

MB 


tfc 

loi 

195 

195 

199 

198 

MS 

Mr 


tr 

105 

ISZ 

196 

199 

IB8 

196 

Me 


ts 

/OS 

ISO 

199 

Mi 

186 

196 

M8 


t 7 

?8 

136 

199 

133 

/ee 

199 

M7 


tio 

IDI 

IS? 

179 

133 

/68 

leo 

/ez 


til 

/O! 

136 

I8B 

131 

170 

166 

186 


tiz 

IDI 

13? 

186 

/39 

I7Z 

/87 

/88 


tt3 

/OS 

197 

190 

MS 

/SO 

199 

MS 


tli^ 

/OS 

196 

/Br 

M7 

/B6 

MS 

Me 


tl5 

79 

63 

78 

89 

78 

7$ 

Bo 


U 

/OO 

191 

199 

MO 

IBB 

/99 

M6 











h-h, 

905 

99 

S9 

?e 

lOZ 

207 

330 



j03 

198 

I99.S 

! 96.5 

188.5 

/?S.S 

MB 


R TUP 8R6 

ess 

77- 

168 

156 

Z99 

723 

/08S 


R eOT 

1010 

— 

Z9Z 

— 

39/ 

9/3 

1307 



2C,40 

9?o 

9?Z 

992 

M9I 

3770 

6060 



1913 

ZDZO 

ZO/0 

Z/90 

Zoso 

20/0 

Zois 





* 

1 



1 







1 


i 

1 

! 



\ 



I 

1 






L - 








_ 



[ 





C.9 












































































DATA SHEET CALCULATED DATA 


RUN NO. 

(o ^ 

&>5 

(o6 

67 

GB 

(>9 

70 


RPn 

8 ZOO 

6200 

6200 

tOj^OO 

f 0,900 

iOj^Too 

\ / 0,900 


p, 

^65 

USO 

l(o/5 

2 oo 

2 05 

9^75 

GB6 



^10 

Ooo 

$50 

f03 

797 

2 73 

i 935 


k =A 

60 

60 

SO 

6/ 

SO 

SO 

99 


FLOW m 

' 16 

26 

26 

10 

76 

26 

! 26 


mfSF 

hZLO 

0-S7.Z 

0 89.8 

^-/z.o 

3-03.5 

2-09,0 

I-/7.7 


rinEffigo 

hOG.Z 

096.7 

0-3 6. Z 

3-0/.0 

2-Z/.0 

1-90.0 

1 hos.g 


t, 

ns 

199 

199 

I9i> 

^90 

795 

198 


^2 

/% 

197 

198 

199 

195 

^ 793 

197 


ts 

ny 

197 

197 

196 

196 

799 

* 197 


t4 

196 

ifd 

I9B 

)95 

196 

196 

197 

1 

ts 

197 

198 

199 

I9L> 

)96 

796 

198 

1 

t(, 

197 

/9e 

197 

196 

195 

799 

197 


tr 

197 

198 

199 

I9(; 

796 

796 

198 


te 

197 

198 

199 

• I9i> 

I9G 

196 

/98 


t? 

I9i> 

195 

196 

l?3 

192 

793 

196 


110 

177 

)8£ 

jBe 

763 

183 

7 8<f 

763 


til 

189 

188 

lee 

186 

/ 65 

192 

7 86 


tiz 

183 

IBB 

189 

J86> 

/ 8^ 

79/ 

786, 


tl3 

I96> 

198 

I9B 

— 

790 

799 

798 


tif 

197 

198 

198 

189 

193 

)99 

/98 


tis 

e(> 

es 

86 

89 

89 

89 

es 


tlb 

19S 

196 

196 

/92 

793 

793 

796 











A-P 

360 

SSO 

Boo 

62 

97 

223 

386 



177 

198 

198 

196 

195.5 

199.6 

197.5 

1 

R TOP BRG 

H30 

i&oo 

2080 

796 

Z96 

723 

U90 


R BOT aflfr 

I3BS 

1882 

2^00 

799 

366 

906 

1336 

1 


6>S8o 

9B50 

I9p50 

968 

1798 

9060 

7ozo 

1 

i 


20/0 

1998 

1985 

2560 

2560 

2590 

2990 






1 




' 
























r ' 



















0 

_° ’ - • • * 




















































































WTA SHEET CALCULATED DATA 



RUN NO. 

7/ 

1 

7^ 




• 

m 


RPtI 

10,400 

10,400 







p, 

1040 

isso 







Pi 

S70 

635 







P =A 

49 

Sf 




1 



FLOW (Lb.) 

ZS 

ZS 




1 

-1 



TIMF 

1 II IL TOP BPt 

1-04.0 

0-49:2 




1 



TIMEffl'iW 

D-S3.0 

0-39.8 







t, 

I9Z 

198 








m 

197 







is 

190 

ns 








I9Z 

nr 





\ 


ts 

I9z 

197 







tt 

l9o 

/94, 







tr 

/9l 

198 







is 

/9Z 

198 







t? 

189 

197 







tio 

/B3 

ns 







in 

187 

188 







tiz 

IB7 

187 







tl3 

I9Z 

t7e 





• 

1 

t;+ 

19 Z 

ne 






1 

tis 

84 

8s 







tif, 

189 

194 





-- 











A-A 

5Zi 

784 






1 

t7tt^ 

19/ 

/96.S 


. 





R TOP 6RG 

i40/ 

/ezz 







R aoT Bflfr 

I6>f8 

ZZ90 








9t/0 

/jj&eo 








ZSOO 

Z4/0 









1 












































7J 

























































































